Polarized Raman study of the phonon dynamics in Pb(Mgi/3Nb2/3)03 crystal 
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Pb(Mgi/3Nb2/3)03 is one of the simplest representatives of the class of lead relaxors and often 
serves as a model system for more complicated compounds. In this paper we analyse both polarized 
and depolarized Raman scattering spectra, measured in the temperature range between 1000 and 
100 K, by multiple peak decomposition. Basing on this analysis, we present the picture of transfor- 
mations in the crystal and the nature of the spectral lines. According to our model, the formation 
of the Fm3m symmetry in the chemically ordered regions as well as the appearance and freezing of 
the polar nanoregions are the consequences of the same phenomenon: ion off-centered displacements 
and their reorientational thermal motion. The lowering temperature progressively limits this motion 
from eight to four, two and only one site. Raman scattering exists due to the presence of the Fm3m 
and polar clusters. It is a product of the multiphonon interaction process, mediated by the presence 
of the dynamic and static disorder. 



I. INTRODUCTION 

Pb(Mgi/3Nb2/3)03 (PMN) is one of the earliest mem- 
bers of the growing family of lead ferroelectric relaxors, 
which, due to its application potential, has been attract- 
ing attention of researchers for many years. However, 
despite significant efforts, many properties of lead com- 
pounds are still unclear and the very nature of the relaxor 
behavior has yet to be understood. The difficulties stem 
from the high complexity of these materials with a high 
degree of compositional, chemical and structural disor- 
der. 

Pb(Mgi/3Nb2/3)03 has a perovskite structure with, 
on average, a cubic Pm3m symmetry, with Mg^+ and 
Nb^+ ions interchanging on B sites. Measurements of 
the dielectric constant show that the crystal does not 
undergo a sharp transition to a ferroelectric phase. In- 
stead, the dielectric constant exhibits a broad maximum 
at Tmax ~ 270 Ki. Below T^ ~ 620 K, it deviates 
from Curie- Weiss law^ and, at T<350 K, exhibits strong 
dispersion^'. The structure of PMN is not homogeneous. 
If Mg and Nb occupation ratio were uniformly and ideally 
equal 1:2 throughout the crystal, it would have a rhombo- 
hedral symmetry R3m. But, a network of superstructure 
clusters with face centered cubic symmetry (FmSm), de- 
stroys the picture. The size of these clusters is of the 
order of 2-3 nm; the distance between the centers of the 
neighboring clusters is near 2.5 nm^. 

The most widely accepted hypothesis, so-called "space- 
charge model" , associates these Fm3m clusters with an- 
tiferroelectrically ordered^ areas where Mg:Nb occupa- 
tion ratio is equal to 1:1, i.e. with local composition 
Pb(Mgi/2Nbi/2)03^. Such 1:1 clusters carry an exces- 
sive negative charge, which is compensated by the posi- 
tively charged host matrix. Network of these charges acts 
as a source of random fields and may prevent the occur- 
rence of the "normal" phase transition. Alternatively, the 
Fm3m clusters might be explained by the local order in 



the form oiVh[Bl+^^\+^]i/2\Q'^+]i/20-^-. This hypoth- 
esis avoids assumption of the considerable space charges 
and calls for other reasons to explain relaxor behavior, 
including the random fields due to the disorder-induced 
local anisotropjiii, effects related to different ionic radiii^ 
or botli2i. 

The presence of areas with special order causes crys- 
talline structure to change in an intricate way. At the 
high temperature limit, the crystal has a uniform primi- 
tive cubic (Pm3m) structure. Lowering the temperature 
leads to distinguishability between Mg^+ and Nb^+ sites 
and to formation in the ordered regions of a face centered 
(Fm3m) superstructure. The temperature of the Fm3m 
clusters formation has yet to be determined. Husson's 
data suggests they presence at T^850 Ki^. Our data 
(see below) raise the upper limit to a temperature higher 
than 1000 K. 

Cooling the sample leads to the nucleation of polar 
nanoregions (PNR's) characterized by the local distor- 
tions. These distortions are responsible for the deviation 
of the dielectric constant from Curie- Weiss law at T(f^. 
Their presence also causes a non-linearity in the tem- 
perature dependence of the optical refraction indesii^. 
Size of the PNR's is equal to 30-80 A depending on 
temperaturai^. The appearance of polar distortions is 
a consequence of the ion off-centering that is a com- 
mon property of many perovskites. By structural sim- 
ilarity, one could expect that the mechanism of the po- 
lar regions formation in PMN is analogous to the one in 
KTai_:rNb^03 (KTN>ii. However, unlike KTN where 
Nb^+ is the only off-centered positive ion, PMN repre- 
sents a more difficult case. Here, not only both ions 
of B-type, Mg2+ and , are shifted by -0.1-0.2 A, 

but also Pb^"*", ion of A-type, is off-centered by ~0.25- 
0.35 Ai2ii2*2aiSL2^. The B-ions tend to shift along a (111) 
direction, thus having eight equivalent positions. The 
shifts of the Pb^+-ions exhibit a spherically symmetric 
distribution^^. It is often assumedii*i2iiii4 that Nh^+ , 
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TABLE I: Interpretations of Raman scattering spectrum from PMN offered by various researdi groups. Note, many of the 



entries in tlie table 


mutuaUy exclude each other. 


Research 
Group 


45 cm"^ 130 cm~^ 260 cm"^ 420 cm"^ 500-600 780 cm"^ 

cm^^ 


HussorsW 


Pb-O stretching modes 0-B-O Mg-O-Mg Nb-O...Nb Nb-O...Mg 

ULliLliilgi & Li ti LUliili^ b Li t; LUiiiitj b Li t; Ll^iiiiigj 


Dimza^ 

(■\QQ')\ 


Coexistence of tetragonal 
aiiLi 1 lioiiiuuiicoi ai piiaoco 


ivianiiKL iryyo) 




ldink*a 

(1994) 


All spectrum is due to one phonon processes, originates 
from rhombohedral symmetry clusters and consists of 9 modes 


Marssi2i^&2aa 

(1996-98) 


TOi TO2 TO3 TO4 TO2+TO4 or LO4 (1996) 

Nb-O-Mg or LO4 (1998) 


Jiang2^ (1999) 


Due to positional disorder of lead Due to Fm3m clusters and vibration of oxygen octahedra 


T 1 :i 33.34 

(1999) 


All spectrum originates from Fm3m clusters, similarly to the case of PbSci/2Tai/203 
VV-Ej?? 

VII-2I23 Alg 

or Disorder-induced scattering from silent modes Breathing mode 
VV and one VH and second-order scattering of oxigen 
have fractal 
character 
another VH-T29 


Dujovne^^ (2002) 


TA+disordcr 



due to its position in the cell and small radius, acts as 
the main ferroelectric agent. 

At high temperature, the ofF-centered ions are free to 
reorient among all allowed ofF-centered positions. Cool- 
ing down and the growing role of the electric interac- 
tions causes the appearance of correlated regions con- 
sisting of several neighboring cells. Each such region is 
characterized by a giant electric dipole moment and a 
local distortion from the cubic symmetry. The PNR's 
are capable of reorientational motion as whole units. In 
KTN, the development of these regions leads, at some 
critical temperature, to a phase transition. However, in 
PMN, the polar regions develop under constraints im- 
posed by the local anisotropy. As a result, no ferroelec- 
tric phase transition, but a glasslike freezing of the polar 
regions is observedii. A bias electric field of magnitude 
~ 1.8 kV/cm, applied in the (111) direction, is sufficient 
to compensate for the influence of the frustrating effects 
and to induce, at T^o ~ 210 K, a transition to the rhom- 
bohedral R3m ferroelectric phase^'^. 

Raman spectra of PMN have been measured and re- 
ported by several research groupsi^i^. But, the high 
complexity of the processes caused difficulties in their 
interpretation, especially due to possible coexistence of 
clusters with Pm3m, Fm3m and R3m symmetry. It has 
been shown that the major spectral lines are of first- 
order character—. In cubic crystals first-order Raman 
scattering is prohibited by symmetry. Explanation of 
the appearance of first-order scattering in principle is 



the one side of the problem. Another side is the as- 
signment of particular lines. The question of the soft 
ferroelectric mode (TOi) and the strong Raman line at 
^--^45 cm~^, which is close to where the TOi is expected 
to be found, has been like a stumbling block for inves- 
tigators. Interpretations of the PMN Raman spectrum 
offered by various research groups are summarized in Ta- 
ble This table reflects the variety of mutually exclud- 
ing opinions ranging from those that explain the light 
scattering by disorder in general (groups of Hussorki^, 
Krainil^i, MarssiSLSSiSSiSa, Dujovne^) to those that con- 
nect it to the presence of some particular type of disorder, 
the one that can lead to the appearance of clusters with 
symmetry allowing first-order scattering, like face cen- 
tered cubic (Jiang^^, Lushnikov^^'^'', SinjiSSi^SiSLSSi^) , 
rhombohedral (Idink^) or tetragonal (DimzE^i). The as- 
signement of particular lines is even more contradictory. 
Recently made first-principle calculations*^ can confirm 
only few of them. We believe, the arguments might be 
resolved by taking into account the possibility of coexis- 
tence of different mechanisms, each playing a role in the 
scattering processes. 

Recently published low branches of phonon dispersion 
curve o^^i'^^1^^1*^ , measured by neutron scattering, must 
be helpful for the phonon assignment of some of the Ra- 
man lines. According to these data, the strong 45 cm~^ 
line has an energy close to the energy of the zone cen- 
ter (ZC) TOi or the zone boundary (ZB) TA phonon. 
Lowering the temperature, TOi zone center phonon ex- 
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hibits softening, at T^T^ K, due to the drastic increase 
of damping, it disappears. The overdamping continues 
down to Tdo ^ 210 K^^ It also causes significant, almost 
sixfold, broadening of the TA phonon brancb^i. It seems 
reasonable to expect that the Raman line at 45 cm~^ 
is also affected by the described phenomenon. However, 
none of the known Raman reports show any significant 
effects that could be attributed to this overdamping. 

The formation of the lower-symmetry clusters in 
the host crystal is accompanied by relaxational and 
reorientational dynamics. This dynamics should be 
reflected in the light scattering spectrum either di- 
rectly or through interaction with the phonon modes. 
Theoretical considerationgi^ show that internal relax- 
ational motion leads to the appearance of the cen- 
tral peak (CP) in the Raman spectrum. Interac- 
tions of this motion with phonons can cause soften- 
ing of them*?.. These predictions have been tested 
on several types of crystals with internal degrees 
of freedom, like alkali-halide-cyanide (K(CN)a;-KBri_2;, 
K(CN)^-KCli_:,)i2iSaaiS compounds. Recently, we 
have analyzed^ the reorientational motion in ferroelec- 
tric KTai_a;Nba;03 crystals. In this paper we discuss the 
role of the relaxational and reorientational dynamics in 
the scattering processes from PMN crystal. 

In the paper presented we offer a complex analysis 
of temperature dependencies of polarized (VV) and de- 
polarized (VH) Raman spectra from a single crystalline 
PMN sample. This analysis is based on multiple peak 
decomposition. An attempt to apply similar kind of 
analysis has earlier been made by Siny et alu^. But, 
it was restricted to the low frequency VV spectra in a 
limited temperature range (with maximum temperature 
~620 K), concentrating on the central peak (CP) only. 
To obtain the complete picture, we consider the temper- 
ature evolution of the whole Raman scattering spectrum 
(up to 1000 cm~^) in the 100-1000 K temperature range. 
We interpret our results in connection with recently pub- 
lished neutron scattering data and with the concept of 
interaction of phonons with internal relaxational motion. 
As a result of this work, we present our view on the pic- 
ture of structural transformations in PMN crystal and on 
the origin of its Raman spectrum. 



II. EXPERIMENTAL SETUP AND RESULTS 

We have investigated the Raman scattering from (100) 
cut PMN single crystal. The crystal was grown by the 
high temperature solution technique using 30 wt% PbO 
as flux. The growth conditions were optimized based on 
the pseudo-binary phase diagram established for PMN 
and PbO'"'^ . The as-grown crystal exhibits a pseudo-cubic 
morphology with (100) cub growth steps. A (100) / (110) 
cub-oriented crystal of a volume of 53 mm'^ was cut from 
a large as-grown crystal and polished with fine diamond 
paste (down to 0.25 mm). It showed very high optical 
quality and satisfied the requirements for light scattering 
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FIG. 1: Examples of Raman scattering spectra measured at 
various temperatures without polarization analysis and cor- 
rected by the population factor. 



studies. The scattering was excited by propagating in 
(100) direction 514.5 nm light from a 200 mW Ar+-ion 
laser, focused to a 0.1 mm spot. The scattered light was 
collected at an angle of 90° with respect to the incident 
beam (i.e., in (010) direction) by a double-grating ISA 
Jobin Yvon spectrometer equipped with a Hamamatsu 
photomultiplier R-649. For most of the measurements, 
the slits were opened to 1.7 cm~^. However, in order to 
acquire more precise data in the central peak region, at 
the temperatures close to the maximum of the dielectric 
constant ( 100 < T < 350 K), the slits were narrowed 
to 0.5 cm~^. Each polarization of the scattered light, 
{x\zz\y) (VV) and {x\zx\y) (VH), was measured sepa- 
rately. In order to exclude differences in sensitivity of the 
monochromator to different polarizations of the light, a 
circular polarizer was used in front of the entrance slit. 
For control purposes, we also took measurements without 
polarization analysis. Finally, to protect the photomul- 
tiplier from the strong Raleigh scattering, the spectral 
region from -4 to -|-4 cm^^ was excluded from the scans. 
The data were collected in the temperature range from 
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FIG. 2; Examples of Raman scatteting spectra measured in VV (a) and VH (b) geometries upon cooling from 1000 to 100 K. 



1000 to 100 K. The cooling rate was 0.5-1 K/min. Every 
50-20 K the temperature was stabilized and the spectrum 
recorded. We agree with KeiM that the reproducibility 
of the spectra in the temperature range approximately 
between 200 and 350 K was not very good. In this range, 
the spectra exhibited a strong dependence on the cooling 
rate and other experimental conditions. 

Figd presents examples of light scattering spectra 
measured at different temperatures without polarization 
analysis. To facilitate comparison, we corrected the spec- 
tra by the Bose population factor: 

F{f T) = I ^'^^ Stokes part 

' 1 for anti-Stokes part ' 

where 

n{f)^{exp{hf/kT)-l)-\ 

Figure [3 demonstrates scattering spectra measured in 
VV (a) and VH (b) geometries at different temperatures. 



These spectra we present in uncorrected, "as-measured", 
form. 

As can be seen from Figs. n]and[21 our spectra are con- 
sistent with those from Refsii^ and^i. In the high tem- 
perature region, a typical spectrum consists of two strong 
lines centered approximately at 45 cm^^ and 780 cm~^ 
(labelled A and B) and of three broad bands (C, D, E). 
Line A exhibits fine structure, which is explicitly seen 
on the spectra measured with polarization analysis (see 
Figs. [5] and 13). Lowering temperature leads to the split- 
ting of the broad bands C, D and E into a number of 
narrower lines (labelled by indices) and to appearance of 
new line F. 

As seen from Fig. ^ starting from 1000 K, the reduced 
intensity of the scattering, first grows until the temper- 
ature of 700 K then sharply decreases to a minimum 
located close to the Burns temperature T^ « 620 K, after 
which it increases again. At ~ 550 K its strength gets re- 
stored. Below this temperature, the intensities of the ma- 
jor lines (especially A and B) are determined primarely 
by the Bose factor. I.e. being corrected by it, they do 
not show significant temperature changes, demonstrat- 



5 




FIG. 3; Examples of multiple peak decomposition of Raman spectra measured at 1000 and 110 K both, in VV and VH geometry. 
For convenience, major peaks are labeled with letters. Inserts show examples of fits for the line A using the model of the two 
coupled harmonic oscillators, as described in Discussion. 



ing characteristic first-order behavior. The shape simi- 
larity of the low- and high- temperature spectra shows 
that even at 1000 K the scattering has first-order char- 
acter. In the perfecly cubic crystal, first-order scattering 
is prohibited. Therefore, one should assume the presence 
of distortions in the form of lower symmetry clusters, 
like Fm3m or R3m. At high temperature, these clusters 
might be present in the dynamic form, having short life- 
time. Lowering the temperature, their lifetime increases 
and they become progressively more static. The observed 
around T^ decrease of Raman intensity is also clearly seen 
in Fig. 12 However, there are no reports on this effect 
in the literature; it might have been overlooked. This 
intensity drop indicates worsening of the optical quality 
of the crystal. In our opinion, it is consistent with the 
formation of Fm3m clusters in their evolution from dy- 
namic to static form. The splitting of peaks is associated 
with the development of R3m clusters. The first signs of 
the splitting appear at T^. Strarting from T/ ~ 350 K 
(which is also close to where the dispersion of the di- 
electric constant appears), this effect becomes even more 
definite and the intensities of the components begin to 
grow. Below Tdo ~ 210 K, which is the temperature of 



electric field-induced ferroelectric phase transition^'^ , the 
shape of the spectra is relatively stable. 

A comprehensive analysis of our data and their com- 
parison with results of other experiments allows to under- 
stand the meaning of the mentioned above temperatures 
Td « 620 K, T/ « 350 K and T<jo ~ 210 K in the tem- 
perature evolution of the Raman spectra and the lattice 
structure. 



III. MULTIPLE-PEAK DECOMPOSITION OF 
SPECTRA PROCEDURE. 

To be able to analyze the data, we decomposed the 
measured spectra using multiple peak fitting procedure. 
It was found that satisfactory fit could be achieved with 
the assumption that the central peak has a Lorentzian 
shape and that each of the other peaks is described by 
the spectral response function, i.e. damped harmonic 
oscillator, modified by the population factor 

(f2jf2fi^2r, mT), (2) 
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where Ti and foi are the damping constant and the mode 
frequency. 

As ah of the peaks are much better resolved at low 
temperatures, we started our fitting procedure at the 
low-temperature end of the data set (at 110 K) and, 
then, observed the evolution of the peaks with increas- 
ing temperature (i.e., in order, opposite to the order of 
measuring). At the same time, we tried to minimize the 
number of peaks necessary to achieve a reasonably good 
fit. The control data set (measured without polarization 
analysis) has been used to calibrate the positions and 
widths of the weak and poorely resolved peaks from the 
VV and VH data sets. Since a large number of parame- 
ters is involved, results of a particular fit may depend on 
their initial values. To stabilize the results, the best-fit 
values of parameters obtained at the previous temper- 
ature were used as initial values for the next one. In 
this manner, several sets of fits were obtained and an- 
alyzed. It is remarkable, that in all of them the major 
parameters showed the same trends of behavior, confirm- 
ing the importance of the mentioned above temperatures: 
Trf 620 ±50 K, T/ w 350 ±25 K and Tdo ~ 210 ±25 K. 
These trends are summarized in Table Hll One of the sets 
have been selected to report the most interesting results 
in detail. Examples of the fits at the temperatures of 
1000 and 110 K in VV and VH geometries are shown in 
Fig. 13 For clarity, we describe the observed phenomena 
from high to low temperatures, following the same order 
as in measurements (unless the opposite stated). 



IV. MOST INTERESTING RESULTS OF 
DECOMPOSITION AND THEIR ANALYSIS. 

A. Central peak. 

The temperature dependences of the fitting parameters 
for the central peak (CP) are presented in Fig. 0] Circles 
correspond to the VV and triangles to the VH component 
of the peak. The existence of the CP is a direct conse- 
quence of the lattice relaxations, which are very sensitive 
to the restrictions imposed by the low-symmetry clusters. 
If relaxations are fasl^, the CP is low-intense and broad, 
whereas their slowing causes growth and narrowing of the 
peak. 

We would like to point out a striking similarity of 
the temperature behavior of the CP in PbMgi/3Nb2/3 03 
(Fig. EI) and in KTao.85Nbo.15O3 (Fig.3 in Ref.^^) gj.yg_ 
tals. We have shownS that the temperature behavior of 
the CP in KTN can be explained by the model involv- 
ing the relaxational motion of off-centered Nb ions and 
its progressive restriction with temperature decrease. In 
the cubic phase, Nb ions are allowed to reorient amongst 
eight equivalent <111> directions. The appearance of 
the PNR's followed by a sequence of phase transitions 
down to a rhombohedral R3m phase, limits the ion mo- 
tion to four, two and, finally, locks it in only one site. 
This model is in agreement with the diffuse neutron scat- 
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FIG. 4: Temperature dependencies of the intensity (a) and 
half-width (b) of the Lorentzian approximation for the cen- 
tral peak in VV (circles) and VH (triangles) geometries of 
experiment. 



tering studies of similar system.^^. The similarity of the 
CP behavior in these two very different materials, sug- 
gests that the temperature evolution of the polar clusters 
in PbMgi/3Nb2/3 03 passes through the similar stages as 
in KTao.85Nbo.15O3 (KTN), however, it is not accompa- 
nied by appearance of the long-range order. 

Starting from the high-temperature end (Fig. 0)), the 
first important feature is the strong and narrow scat- 
tering in the VV geometry accompanied by a relatively 
weak scattering in the VH geometry. This indicates the 
presence of a symmetric slow relaxational motion, involv- 
ing 180° reorientations of ions. Lowering the temper- 
ature, starting from ^ 900 K, a cessation of this mo- 
tion causes an intensity decrease of the CP. Slow at first, 
this decrease becomes sharper and reaches minimum near 
Trf « 620 K. The prohibition of I8O0 reorientations in- 
troduces anisotropy into the lattice causing the distin- 
guishability between Mg and Nb occupied cites and to 
formation in the 1:1 ordered areas of the superstructure 
with, on average, Fm3m symmetry. It also imposes first 
restrictions on the reorientational motion of the dynami- 
cal R3m polar nanoregions. Now, they can reorient only 
amongst four neighboring <111> directions, forming, on 
average, tetragonal-like distortions. These processes are 
accompanied by the appearance of large (of the order of 
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TABLE II: Major peaks and their response to the temperature restrictions on ion motion. Only changing properties are shown. 
Td ~ 620 K is also characterized by the total loss of intensity of scattered light. 
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VV 


Hardens and narrows 


Fine structure appears 


Splits in two 




VH 
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VH 
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Peak F 
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Peak E 
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Fine structure appears 


Splits in two, VH appears 


Peak B 


VV 


Hardening 


Hardening 


VH appears 



wavelength of light) dynamic fluctuations causing wors- 
ening of the optical quality of the sample. 

With further decrease of the temperature, the opti- 
cal quality of the crystal improves. From ~ 550 K, the 
four-cite reorientational motion of the PNR's starts to 
slow down, which is marked by the narrowing of the 
VV component of the CP and increase of its VV and 
VH intensities. The simultaneous broadening of the VH 
component up to the maximum at ~ 450 K, indicates re- 
arrangements in the crystalline structure leading to the 
appearance of the new type of restrictions on the ion mo- 
tion. Analogy with KTN suggests that below ^ 450 K 
the motion of R3m clusters becomes limited to two neig- 
boring <111> orientations, averaging in monoclinic-like 
distortions. Such a rearrangement causes some decrease 
in intensity of the VV component (with minimum at 
^ 400 K), while the VH intensity keeps growing. Be- 
low ^ 400 K, the slowing down of the two-site relax- 
ational motion causes narrowing the CP and increase of 
intensities of both components. From T/ « 350 K, these 
effects become especially dramatic. Further temperature 
decrease, starting from ^ 300 K, leads to the complete 
prohibition of intersite reorientational motion of R3m 
clusters, i.e. to appearance of static R3m clusters. This 
is marked by a sharp decrease in intensity of both com- 
ponents of the CP. At the temperature T^o ~ 210 K, the 
process of conversion of PNR's from dynamic to static 
form (freezing) is primarily finished. Below this temper- 
ature, the central peak is narrow and its intensity is small 
in both scattering geometries. 

We should mention that the result of our analysis of 
the VV component of the CP exhibits similar tendencies 
to those in Refi^, except at the low temperature end. 



The discrepancy occurs due to the difference in experi- 
mental technique. In the central peak region, we did the 
measurements with much higher resolution (0.5 cm~^, as 
compared to 2 cm~^ from RefjM). Consequently, we were 
able to provide a more correct separation of the slow re- 
laxations from the elastic scattering, which was especially 
important at low temperatures. 



B. Line A. 

Figure |51 presents the temperature evolution of the fit- 
ting parameters of the peak A, showing its position (a), 
reduced intensity (b) and damping constant (c). This 
peak has a triplet structure, containing one component 
in VV (circles) and two components in VH (up and down 
triangles) geometries. The fitting parameters for this 
peak exhibit changes at mentioned above temperatures 
Trf, T^ and T^o, confirming their importance for the 
structural evolution of the crystal. However, the origin 
of this peak (see Table. P) requires qlarification. From 
a comparison with the frequencies of the phonon modes 
determined from neutron scattering, it is clear that this 
line cannot be due to the zone center soft TOi mode 
(black stars in Fig. |3J). On the other hand, the lower 
frequency VH component, and possibly the VV compo- 
nent could be due to the disorder-induced zone boundary 
scattering on TA phonon (white stars in Fig. |SJ|. How- 
ever, the higher frequency VH component would still not 
be accounted for. 

In an attempt to account for both VH components 
simultaneously, we have tried to make use of a coupled 
oscillator model. In the approximation of linear coupling 
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FIG. 5: Temperature dependences of the fitting parameters: 
position (a), reduced intensity (b) and damping constant (c) 
for the triplet line A (located at around 45 cm~^). Phonon 
frequencies, measured by neutron spectroscopy— , are 
shown for comparizon. 

between two harmonic oscillators^^ this approach gains: 

^ r / N 11 fCY-BZ\ 
liuj) = [niu;) + 1] ^ j^2 + c^ ) ' 

where 

B = (ujf -UJ^){UJ^ ~UJ^) +UJ^{jf2 -7172) - WlW2Af2, 
C = Lu[-fl{u;^ - LU^) + 72(^1 - UJ^) ~ 2(wiW2)^/^Al27l2], 
Y = 5'2(c^2 _ ^2) ^ ^2(^2 _ ^2) _ 2^i^2K0^2)^/'Ai2, 

Z = u;iSh2 + Shi - 2^1 52712). 

In these equations lui and 0-12 are the resonant frequencies, 
71 and 72 are the damping constants. Si and ^2 are the 
structure factors for two oscillators, A12 and 712 are the 
real and imaginary parts of the coupling constant. 

It was shown^, that a unitary transformation can be 
found, which would set either real or imaginary part of 
the coupling to be equal to zero. We have chosen to 
use imaginary coupling and set A12 = 0. We also set 
the structure factors and S2 to be the same at all 
temperatures. The rest of the parameters were treated 
as being temperature dependent. In the result of fitting, 
we have found that the values of structure factors can be 



FIG. 6: Temperature dependences of the fitted parameterst 
of the line A VH component (located at 45 cm~^) using the 
model of coupled harmonic oscillators. From top to bot- 
tom: resonant frequencies (a), damping constants (b) and 
coupling constant (c). Phonon frequencies, measured by neu- 
tron spectroscopy^^*— are shown for comparizon. 

approximated as 5i = 22.3 cm~^ and S2 = 20.6 cm~^. 
The other parameters as a function of temperature are 
shown in Fig|BJ the best-fit frequencies (a) , the damping 
constants (b) and the coupling coefficient (c). 

The examples of the fits (inserts in FigO demon- 
strate that this model gives a good approximation to 
the shape of the peak A. However, a comparizon of the 
best-fit frequencies (Fig. Eli) with the neutron scatter- 
ing datai^Si^^i^^*^ rules out the possibility to explain the 
higher-frequency VH component by the interaction be- 
tween ZC TOi (black stars) and ZB TA (white stars) 
modes. The lower frequency oscillator can still be asso- 
ciated with disorder-induced scattering on ZB TA mode. 

The possible explanation of the presence of two VH 
components in the peak A can be in the correspondence 
to TA modes with different polarizations in different crys- 
talline planes of a distorted lattice. (A similar effect we 
have seen in monodomain KTai„a;Lia;03, where one or 
the other TA polarization was observed, depending on 
the orientation of the electric field.) These differently 
polarized modes are influenced by the polar nanoregions, 
which explains the increase of damping (Fig. Eb) and 
coupling (Fig. Et) coefficients with lowering tempera- 
ture. 
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It is also necessary to point out an observation con- 
cerning the temperature dependences of the frequencies 
of the components of the triplet A in approximation of in- 
dependent oscillators (Fig. [S^). These dependences seem 
to show anomalies at intersections with ZC TOi modes 
(black stars). For example, the VV component exhibits 
deeps at ~ 700 and 350 K and VH-I component has a 
knee at ~ 400 K. 



13 It also exhibits similar temperature tendency to that 
described in Ref.^^. Detailed understanding of this phe- 
nomenon requires further work, and it is a subject of 
separate paper. 
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FIG. 7: Amplitude of splitting between the components of 
the line E. The insert demonstrates that this splitting follows 
the critical dependence. 



C. The most important parameters of other lines. 

Among the changes occuring to the Raman spectrum 
with cooling the sample, the splitting of the phonon lines 
desrves special attention, since such a splitting is a di- 
rect indicator of the modifications of the local structure. 
As we already mentioned, the splitting occurs to the 
lines C, D, and E, which is especially well seen in the 
VV geometry (Fig. At high temperatures, their 

widths are very big, so the lines almost merge in a sin- 
gle broad shoulder. Lowering the temperature, they be- 
come more discernible and the fine structure gradually 
develops. The first signs of it appear at T^ « 620 K. 
Below Tpr « 350 K, the fine structure becomes more ob- 
vious. These observations support the expressed above 
idea that the scattering originates from the distorted re- 
gions. Highly dynamic and disordered at high tempera- 
ture, they become progressively more static with cooling 
down, and their motion becomes more correlated. This 
process is reflected by the temperature evolution of the 
lines C, D, and E. 

Below the temperature of Tpr ~ 350 K, the splitting 
value of the line E exhibits especially interesting trend. 
As seen from Fig. \7\ with temperature decrease be- 
low Tpr the distance between the E-line components in- 
creases. Moreover, this increase exhibits the characteris- 
tic order parameter behavior, following the (Tp^— T)~^/^ 
dependence, which is demonstrated in the insert to Fig. 
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FIG. 8: Temperature dependencies of the fitting parameters: 
position (a), reduced intensity (b) and damping constant (c) 
of the line B (located at ~780 cm^^.) 



Finally, we would like to stop on the line B located at 
the upper end of the spectrum. Fig. |S1 demonstrates the 
fitting paramenters for this line: position (a), reduced in- 
tensity (b), and damping coefficient (c). This line is very 
strong in the VV geometry but it is rather weak in the VH 
one (i.e. characterized by Aig symmetry). It seems to be 
the least influenced by the ordering processes in the crys- 
tal. However, the temperature evolution of its parame- 
ters reflects development of the low-temperature phase, 
confirming the expressed above ideas about formation 
of dynamic polar clusters at T^ « 620 K, their slowing 
down, and appearance of the static order at Tp^ ~ 350 K. 
As it is seen from the Table n the appearance of this line 
in the Raman spectrum has several contradictory expla- 
nations. However, according to the recent first-principle 
calculations^, this mode most likely is a consequence of 
oscillations of oxygen ions. 
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V. CONCLUSION 

This paper offers a comprehensive analysis of the Ra- 
man scattering spectra of PbMgi/3Nb2/3 03 in the tem- 
perature range between 100 and 1000 K. From our analy- 
sis, the structural evolution of the lattice occurs in several 
stages. Even at 1000 K, the PMN crystal is characterized 
by the presence of dynamic distortions of the lattice from 
cubic symmetry, which are responsible for existence of 
first-order Raman scattering at such a high temperature. 
At the Burns temperature « 620 K, the first restric- 
tions on the rcorientational motion of the off-centered 
ions appear. These restrictions lead to slowing of the 
motion and growth of correlations between neighboring 
regions with R3m symmetry. They also introduce distin- 
guishability between sites occupied by Mg and Nb, which 
causes onset of Fm3m symmetry in the 1:1 ordered areas. 
Further cooling leads to progressive growth of restric- 



tions on ionic motion and, starting from Tpr ~ 350 K, 
to the appearance of the static polar nanorcgions. At 
Tdo ~ 210 K, the formation of the low-temperature phase 
is primarily finished. The shape of the phonon lines is 
determined by multiphonon processes involving phonons 
with different polarizations propagating in different di- 
rections, which interact with dynamic and static disor- 
der. 
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